Anisotropic pair-superfluidity of trapped two-component Bose gases 
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We theoretically investigate the pair-superfluid phase of two-component ultracold gases with neg- 
ative inter-species interactions in an optical lattice. We establish the phase diagram for filling n = 1 
at zero and finite temperature, by applying Bosonic Dynamical Mean-Field Theory, and confirm the 
stability of pair-superfluidity for asymmetric hopping of the two species. While the pair superfluid 
is found to be robust in the presence of a harmonic trap, we observe that it is destroyed already by 
a small population imbalance of the two species. 

PACS numbers: 03.75.Kk, 03.75.Mn, 67.85.-d, 67.85.Hj 
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I. INTRODUCTION 

Ultracold gases in an optical lattice are promising 
as flexible quantum simulators for the study of quan- 
tum phases which are not easily accessible in condensed- 
matter physics [l|. Experimental realization of the 
supcrfluid-to-Mott transition has paved the way to stud- 
ies of strongly interacting Bose or Fermi gases in opti- 
cal lattices J2J. Recently, Bose-Bose mixtures of 87 Rb 
and K have been produced and loaded into an optical 
lattice [3|, which provide a new playground for investi- 
gating the interplay between kinetic energy, interaction 
and the spin degree of freedom. In a further line of in- 
vestigation, a Bose-Bose mixture of two hypcrfine states 
of 87 Rb has been used as a spin-gradient thermometer 
for measuring the temperature of ultracold gases in an 
optical lattice [J, Q. Moreover, the effect of a second 
species on bosonic superfluidity in an optical lattice has 
been studied [y] . In these mixtures of different species or 
different hyperfinc states, the most fundamental physics 
associated with quantum magnetism and the spin degree 
of freedom can be explored [7[ • 

Recently, it was found that a three-body hard-core 
constraint can stabilize a system of single-component 
bosons with attractive two-body interactions in an opti- 
cal lattice [8[, and numerical simulations have been per- 
formed to study properties of this system [9l-[llj|. For 
two-component bosons with attractive inter-species in- 
teractions in an optical lattice, one interesting ground 
state is the pair-superfluid phase (PSF), which has been 
predicted and studied theoretically in several papers [12| — 
[l9j . Qualitatively, the PSF can be viewed as a con- 
densate of bosonic pairs of different species or hyperfinc 
states due to a second-order hopping of the pairs, but 
with a strongly suppressed first-order tunneling of single 
atoms. For hard-core bosons with total filling n = 1, the 
PSF for attractive interactions is equivalent to the XY- 
ferromagnetic phase (or super-counter-fluid state, de- 
noted as SCF) for repulsive interactions, i.e. the PSF 
consists of particle-particle pairs of different species while 
the XY-ferromagnetic phase is composed of particle-hole 
pairs. It is apparent that the PSF can only exist at very 
low temperature compared to the critical temperature of 



quantum magnetic phases. The latter are also governed 
by second-order tunneling processes, leading to an effec- 
tive spin-exchange coupling [2(J, [21|, which has already 
been observed for a double- well system [22|, [23| . At the 
current stage, it is experimentally challenging to reach 
the critical temperatures of quantum magnetic phases 
and the PSF [f|. It is expected that these phases can 
be detected via momentum-space correlations observed 
in time-of-flight measurements [171 . l24f or by optical mi- 
croscopy with single-site resolution [2a, |26[ . 

Previous studies of the PSF in a two- or three- 
dimensional optical lattice mostly focus on symmet- 
ric parameters for the two species (for an exception 
in [l2|, UM) , since this is the most favorable condition 
for the PSF [lT]]- However, there is still a lack of de- 
tailed quantitative investigations of the PSF for the ho- 
mogeneous system with asymmetric hopping amplitudes 
of the two species, for the trapped system, and for im- 
balanced mixtures of the two species. Here we inves- 
tigate the properties of the PSF of two-component ul- 
tracold bosons with attractive interspecies interaction, 
both in a homogeneous and a trapped optical lattice. 
This system can for sufficiently low filling be well de- 
scribed by a single-band Bosc-Hubbard model with pure 
onsite interaction. We investigate the homogeneous sys- 
tem by means of bosonic dynamical mean field theory 
(BDMFT), which is a non-perturbative approach towards 
strongly-correlated bosonic systems [27| , and the trapped 
system by re al-sp ace bosonic dynamical mean field theory 
(RBDMFT) [28[ , which includes arbitrary inhomogeneity 
such as a harmonic trap. For the homogeneous system, 
we focus on the phase diagram with filling number n = 1 . 
In particular, we also present the phase diagram for the 
experimentally realized heteronuclear 87 Rb - 41 K mix- 
tures, where double-species Bose-Einstein condensates 
with negative inter-species interactions have been ob- 
served [29( . For the trapped Bose-Bose mixture we study 
the coexistence of Mott insulator, superfluid and PSF. 

The paper is organized as follows: in section II we give 
a detailed description of the Bose-Hubbard model and 
the BDMFT approach. Section III covers our results on 
the homogeneous Bose-Bose mixture in an optical lattice 
and the effect of the trap. Wc conclude in Section V. 



II. MODEL AND METHOD 



III. RESULTS 



We consider a two-component bosonic mixture loaded 
into a 2D or 3D optical lattice. In experiments this mix- 
ture could consist of two different species, e.g 87 Rb and 
41 K as in [3| or two different hypcrfine states of a sin- 
gle isotope, e.g 87 Rb as in [J]. Besides the optical lat- 
tice, we also impose an external harmonic trapping po- 
tential which introduces inhomogeneity in the system. 
For sufficiently low filling, the whole system can be ef- 
fectively described by a two-component inhomogeneous 
Bose-Hubbard model within the single-band approxima- 
tion 

H= - y^ tu(b iv bj U + h.c) + - y^ U\ u ni,\(hj\ - 5\ v ) 

(1) 
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where (i,j) denotes summation over nearest neighbor 
sites i,j, and the two bosonic species are labeled by 
the index A(^) = b, d. b\ v (bi V ) denotes the bosonic 
creation (annihilation) operator for species v at site i 
and fii^ v = b iv bi V the corresponding local density. Due 
to possibly different masses or a spin-dependent optical 
lattice, these two species generally hop with non-equal 
amplitudes tb and td- V\ v denotes the inter- and intra- 
species interactions, which can be tuned via Feshbach 
resonances or spin-dependent lattices [0, [30] ■ l±v is the 
global chemical potential for the two bosonic species, and 
Vi = Vo r i denotes the external harmonic trap, where Vq 
is the strength of the harmonic trap and r^ is the distance 
from the trap center. 

To investigate properties of the system, we apply 
BDMFT for a homogeneous lattice |27| , and its real-s pac e 
generalization (RBDMFT) for the trapped system [281 ] . 
Within RBDMFT, the Hamiltonian ((T|) is mapped onto 
a set of individual single-site problems, and the physical 
properties on each lattice site are determined by a local 
effective action which is captured by an effective Ander- 
son impurity model [271 ] solved by exact diagonalization 
(ED) [ll|3l|. RBDMFT fully captures the inhomogene- 
ity of the system and is capable of providing an accurate 
and non-perturbative description of quantum phases and 
their excitations. The detailed formalism is presented 
in HH- 

To include the effect of spatial inhomogeneity, we also 
employ an LDA approximation combined with single- 
site BDMFT. The advantage of this approach is the 
larger system size accessible. Within LDA+BDMFT, 
the local chemical potential for each species is set to 
Hv(r) — ^v — V(r). In this work, we apply both RB- 
DMFT and LDA+BDMFT to the 2D trapped square 
lattice, and only LDA+BDMFT to the 3D trapped cubic 
lattice. 



In this section, we will investigate properties of two- 
component bosonic mixtures with negative inter-specics 
interactions in a homogeneous 3D optical lattice at zero 
and finite temperature, and also in harmonically trapped, 
inhomogeneous systems both in 2D and 3D. For the ho- 
mogeneous system, we will study the stability of the PSF 
against asymmetric hopping and finite temperature. In 
the presence of an external harmonic trap, we will inves- 
tigate finite size effects of the PSF both in 2D and 3D 
within RBDMFT and LDA+BDMFT. We choose the ab- 
solute value of the inter-species interaction U = Ubd as 
our unit of energy in the following. 
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FIG. 1: Zero-temperature phase diagram for two-component 
bosons with attractive inter-species interaction in a 3D cubic 
lattice. The total filling is n = 1 with rib = n<i = 0.5. 



A. Bose-Bose mixture on a 3D cubic lattice 

We first investigate Bose-Bose mixtures with negative 
inter-species interaction U < in a 3D cubic optical 
lattice. The system is unstable for \U\ > Ub,d, since 
the strongly attractive inter-species interaction cannot be 
compensated by repulsive intra-species interactions, lead- 
ing to a collapse of the system. Here we will demonstrate 
the stability of the PSF in the regime \U\ < Ub,d- We 
explore the zero-temperature phase diagram with asym- 
metric hopping parameters at total filling n = 1 (with 
lib = nd = 0.5) for different interactions, as shown in 
Fig. [TJ We observe two different phases: the PSF with 
(b) = (d) = but 4>bd = {bd) - (b) (d) + 1 0, and the su- 
pcrfluid phase (SF) with (b), (d) > 0. In particular, we 
confirm the existence of a PSF for asymmetric hopping 
amplitudes tb + 1 td- In the regime of weak hopping for 
both species, first-order tunneling is suppressed by the 
strong interactions. Formation of bosonic pairs between 
different species can thus be energetically favored and 
will typically compete with singlc-spccies condensation, 
since the bosonic pairs can hop via second-order tunnel- 
ing. As a result, the PSF will have a non- vanishing order 
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FIG. 2: Upper: Finite-temperature phase diagram for two- 
component bosons in a 3D cubic lattice (T = 0.001|£/|). The 
interactions are set to Ut = Vd = 12|t/| and the total filling is 
n — 1 with rib — Ud — 0-5. Lower: Local density-density cor- 
relator g( 2 ' as a function of temperature for different hopping 
values tb = td- Inset: Zoom of main figure for t v — 0.04|[/| 
(y — b,d), where the red dashed line shows the disappearance 
of the PSF correlator <f>bd- 



parameter <p bd but vanishing superfluid order (b) and (d) . 
On the other hand, when both species acquire large hop- 
ping, a superfluid phase will appear with (6) > and 
(d) > 0. Note that, when the intra-species interaction 
U bid decreases from U b = U d = 12\U\ to 2\U\, the PSF 
shrinks due to the decrease of the effective pair-tunneling 
amplitudes. We remark here that the transition from a 
PSF to a superfluid phase for symmetric hopping ampli- 
tudes occurs at the same value of t v /U b d as that from 
the XY-ferromagnetic to the superfluid phase in the cor- 
responding system with U bd > of equal magnitude. In 
contrast to the XY-ferromagnetic phase, however, the 
PSF also exists for non-integer filling. 

The effect of finite temperature is shown in Fig. [2] 
Generally, the PSF is sensitive to temperature, since 
the pairs are formed in the weak hopping regime and 
their coherence can be easily destroyed by thermal fluc- 
tuations, due to their small effective tunneling of order 
0(t 2 /U). At finite temperature, the PSF regime shrinks 
in the weak hopping regime in favor of developing a 
new unordered phase with vanishing values for <p bd , (b) 



and (d). To further understand this unordered phase, 
we calculate the dependence of the local density-density 
correlator g^ 2 > = (n b rid) — {n b ){nd) on temperature, as 
shown in the lower panel of Fig. [2] We observe that 
g' 2 ' starts to decrease noticeably only above tempera- 
tures of the order of 10 _1 |t/|, which indicates that local 
pairs still exist below this temperature. We therefore 
conclude that the unordered phase, shown in the upper 
panel of Fig. [2] consists of non-coherent pairs of differ- 
ent species. In sufficiently deep optical lattices, these 
pairs are localized, while for larger hopping the pairs de- 
localize over the whole lattice. As a result, the local 
density-density correlator decreases as a function of t b ,d, 
as shown in the lower panel of Fig. [2J Another interest- 
ing feature of the temperature dependence of g^> is the 
increasing (non-monotonic) trend at low temperatures, 
since thermal fluctuations first localize and then break 
the pairs. We remark here that the temperature regime 
of non-condensed pairs (« 0.1|J7|) is experimentally ac- 
cessible j32(, and could be detected via radio frequency 
spectroscopy [3J]. Note that the unordered phase ob- 
served here is qualitatively different from that melted 
from the XY-fcrromagnctic phase, since it can also exist 
for non-integer filling. We also observe that the (single- 
particle) superfluid phase remains almost unchanged for 
the temperature considered here. 
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FIG. 3: Critical temperature of the PSF as a function of 
hopping amplitudes % = td on a 3D cubic lattice with total 
filling n = l. Inset: melting of the PSF vs. temperature along 
the vertical dashed line with hopping amplitudes tb = td = 
0.04 \U\. 

One crucial question regarding the observation of the 
PSF is how fragile it is against finitc-tempcrature effects. 
To address this issue, Fig. [3] shows T c as a function of 
the hopping amplitudes t b = td at different interactions. 
We notice that T c rises as the hopping amplitudes in- 
crease, due to the growing second-order tunneling which 
stabilizes long-range order. The inset of Fig. [3] shows the 
temperature dependence of <ft b d- It indicates a second- 
order phase transition from the PSF to the unordered 
phase. We also observe that the critical temperatures for 
the PSF shown here are comparable to those of the XY- 



ferromagnetic phase [28[ and notably smaller than the 
coldest temperatures which have been measured in most 
experiments until now, with the exception of the MIT 
group where temperatures as low as 350pK (ps O.OlUbd 
with tb/Ubd ~ 0.029) have been achieved @. 
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FIG. 4: Upper: Comparison of the zero-temperature phase 
diagram for two-component hard-core bosons with the one 
obtained by a tensor-product-state approximation [19| for 
a square lattice with symmetric parameters: t — tt — £<i, 
fx — m, — fj,d- The red solid lines are phase boundaries ob- 
tained by BDMFT for U b = U d = 200|[/|, while the blue 
dashed lines are the results of the tensor-product-state ap- 
proximation. Lower: Zero-temperature phase diagram for 
two-component soft-core bosons with Ub = Ud = 2|[/| ob- 
tained via BDMFT. 

To verify the validity of the BDMFT results, compari- 
son has been made with a hard-core boson model solved 
by a tensor-product-state approximation [19| . as shown 
in Fig. 2J We find excellent agreement between the two 
methods. We also plot the phase diagram for soft-core 
bosons (Ub = Ud = 2|C/|), and observe that in this case 
the phase boundary between the PSF and the supcrfluid 
phase is shifted to lower hopping values. 



B. Rubidium-potassium mixture 
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FIG. 5: Phase diagram for a mixture of 87 Rb and 41 K in a 
3D cubic lattice as a function of lattice depth s and Rb-K 
scattering length. The total filling is n = 1 with rib — n d = 
0.5. 



for mixtures of hyperfine states of 87 Rb [J]. However, 
this symmetry is not present for mixtures of Rb and 
40 K, where a negative inter-species interaction has been 
achieved via a Feshbach resonance [29| ■ Here we consider 
a mixture of 87 Rb and 41 K loaded into a 3D cubic lattice 
with wavelength A = 757 nm, which yields equal dimen- 
sionless lattice strength s for the two species. Due to the 
different masses, the ratio of the intra-species interaction 
strengths is then fixed to Ur^/Uk = WKflRb/wRbaK ~ 
0.72 and the ratio of the hopping amplitudes to tjib/tK ~ 
0.47, where £r is the recoil energy. 

Now we explore the phase diagram of 87 Rb and 41 K 
mixtures in a 3D cubic lattice and make predictions for 
ongoing experiments. Since the depth s of the optical lat- 
tice and the inter-species scattering length a^bK are tun- 
able with high accuracy, we show in Fig. [5] the phase dia- 
gram in the a^bK-s plane for total filling n = 1 (rib = rid) 
at zero and finite temperatures. At zero temperature, 
two phases appear: superfluid and PSF. When the scat- 
tering length ORbK is small, the system is in a superfluid 
phase for a shallow lattice. When the depth of the lattice 
is increased, the ratio of tRb/^Rb decreases, resulting in 
a strong suppression of first-order tunneling. The domi- 
nant process will then be hopping of composite pairs. At 
finite temperature, the PSF can be easily destroyed by 
thermal fluctuations which induce a second-order phase 
transition into a unordered phase with (f>bd = 0. Since 
the PSF is formed via second-order tunneling and the 
corresponding energy scale is very small, for the param- 
eters chosen here, this transition already occurs at a low 
temperature T = 0.0005.Er. On the other hand we ob- 
serve that at the same temperature the superfluid phase 
is very stable and the phase boundary between superfluid 
and PSF is almost unchanged. 



Our investigations have so far focused on symmetric 
interactions Ub — Ud, which is a good approximation 
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FIG. 6: Density distributions rib (black line), order parame- 
ters 4>b (green line) and PSF correlator 4>bd (red line) vs. radial 
distance r for different hopping amplitudes at zero tempera- 
ture in a 2D square lattice, obtained by RBDMFT. The inter- 
actions are Ub — Ud = 2|C/|, hopping amplitudes t — tb = td 
and harmonic trap Vo = 0.0002|C/|. 



C. Trapped Bose-Bose mixtures in 2D and 3D 
lattices 

In this section, we simulate the two-component bosonic 
system in both 2D and 3D in the presence of a harmonic 
trap, as is relevant for most experiments. In particular, 
we investigate the stability of the PSF in the trapped 
system. Here we choose a 41 x 41 square lattice for the 
2D case and a 41 x 41 x 41 cubic lattice for the 3D case. 
In 2D, in our simulations we apply both RBDMFT and 
BDMFT+LDA, while in 3D we only use BDMFT+LDA 
due to its lower computational effort. 



1. Balanced mixture 

Fig. [S] shows the density distributions n&, order pa- 
rameter fa and correlator (pbd for the PSF vs. radius 
r at different hopping amplitudes in a trapped 2D op- 
tical lattice. Since the PSF is stabilized only within a 
narrow region for the symmetric parameters (see Fig. 2]), 
the harmonic trap should be very shallow and the hop- 
ping amplitudes need to be fine-tuned. Otherwise, the 
system will go through a phase transition directly from a 
Mott-insulating to a superfluid phase. Here we choose 
completely symmetric parameters: t — tb = td and 
Ub = Ud = 2|f7| with balanced filling for the two compo- 
nents. Therefore only one value for nb,d an d 4>b,d respec- 
tively is shown in Fig. [BJ We observe that a wedding- 
cake structure appears in the trapped system, and the 
coexistence of different phases sensitively depends on the 
hopping amplitudes. At lower hopping t = 0.55|C/|, only 
two phases appear, and the corresponding phase transi- 
tion is from a Mott insulator with total filling n = 2 to 
a PSF with total filling < n < 2 indicated by the non- 
vanishing value of (f>bd while <f>b = 0. If we increase the 
hopping amplitudes, the first-order tunneling of single 
atoms will increase, which induces large density fluctu- 



ations in the system, leading to a phase transition from 
the PSF to the superfluid phase. We clearly observe this 
point from the middle panel of Fig. [5] where the superfluid 
phase starts to appear right in the middle of the PSF. 
With further increase of the tunneling amplitudes the su- 
perfluid dominates, as shown in the lower panel of Fig. [51 
where the PSF completely disappears at t = 0.7|[/|. 




FIG. 7: Comparison between results from RBDMFT (R) and 
those from BDMFT+LDA (L) for a 2D square lattice. Den- 
sity distributions rib, order parameters <j>b and PSF correla- 
tor 4>bd vs. radial distance r at zero temperature in a 2D 
square lattice. The interactions are Ub = Ud = 2|[/| and the 
hopping amplitudes tb — td — 0.09\U\ with a harmonic trap 
Vo = 0.0002|[/|. 
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FIG. 8: Density distribution (black line), order parameter 
(green line) and PSF correlator (j>bd (red line) vs. radius r 
for different hopping amplitudes at zero temperature for a 
trapped 3D cubic lattice obtained within BDMFT+LDA. The 
interactions are Ub = Ud — 12|f7|, hopping amplitudes t = 
tb = td and harmonic trap Vo = 0.0002|[/|. 

Fig. [7] shows a comparison between the results of RB- 
DMFT and those of BDMFT+LDA for a 2D square lat- 
tice. We observe good agreement between the two meth- 
ods except close to the phase transition. In spite of the 
artificially sharp phase transition feature of LDA, the re- 
sults of BDMFT+LDA are still reliable with sufficient 
accuracy in the regime away from the transition. We 
will therefore apply BDMFT+LDA to tackle the 3D case 
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FIG. 9: Temperature dependence of the PSF correlator 
as a function of radius r for a 3D cubic lattice obtained 
within BDMFT+LDA. The interactions are U b = U d = 12|C/|, 
hopping amplitudes tb — ta — 0.045|£/| and harmonic trap 
V = 0.0002|t/|. 



due to its higher computational efficiency compared to 
RBDMFT. 

Let us now investigate the stability of PSF in the 3D 
case in the presence of a harmonic trap. Results ob- 
tained within LDA are shown in Fig. 8. Here we choose 
completely symmetric parameters: t = tb = td and 
Ub = Ud = 12|C/| with balanced filling for the two com- 
ponents. And only one value for rib. d and <pb. d is shown in 
Fig. [SJ respectively. Compared to 2D, we observe a simi- 
lar scenario of phase coexistence in the trapped 3D cubic 
lattice: at lower hopping amplitudes the Mott insulator 
and the PSF are coexisting, at intermediate hopping am- 
plitudes, the superfluid phase appears due to increased 
density fluctuations, and at even larger hopping ampli- 
tudes, the PSF will disappear. 

We are also interested in the effects of temperature on 
the PSF. Fig. [9] shows the correlator <j>bd for different 
temperatures. We observe that the PSF is sensitive to 
thermal fluctuations. At finite T, the PSF is reduced in 
favor of developing an unordered phase characterized by 
4>bd = 0. On the other hand, the density distribution is 
rather insensitive to small finite T. 



2. Imbalanced mixture 

As pointed out above, asymmetry in the hopping of 
the two species does not destroy the PSF. However, fill- 
ing imbalance hinders the formation of the pairs [171 . |34| . 
We will now study this effect in more detail. The imbal- 
ance, Nb ^ Nd, will be controlled by a nonzero chemical 
potential difference A/i = /j/, — [id which can be viewed 
as an effective magnetic field. Results in 2D are shown 
in Fig. [TUJ Upon increasing A/i, the PSF will cease to 
exist and be replaced by a superfluid phase, since the 
chemical potential difference eventually exceeds the pair- 
ing gap for the PSF, allowing unpaired excess atoms to 
enter the PSF region. We find that even small population 
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FIG. 10: Density distributions rib,d, order parameters (b), (d) 
and PSF correlator <f) b d vs. radial distance r for different A/i 
at zero temperature in a trapped 2D cubic lattice, obtained 
within BDMFT+LDA. Panel (d) shows the filling difference 
(rib — nd) vs. radius r. The interactions are Ub — Ud — 12|C/| , 
hopping amplitudes tb — td — 0.05\U\, (/it + /id)/2 = 0.48 and 
harmonic trap V = 0.00015|£/|. 




FIG. If: Density distributions n^, order parameters (6), (d) 
and PSF correlator <f>bd vs. radius r for different A/i at zero 
temperature in a trapped 3D cubic lattice, obtained within 
BDMFT+LDA. Panel (d) shows the filling imbalance (n b - 
rid) vs. radius r. The interactions are Ub — Ud — 12|f/|, 
hopping amplitudes tb — td — 0.04|£/|, (fj,b + A*d)/2 = 0.47 
and harmonic trap Vo — 0.0003|f7|. 



imbalance can destroy the PSF in a trapped system. On 
the other hand, as shown in Fig.[T0l the density distribu- 
tion is almost unchanged for similar parameters. When 
increasing the imbalance parameter A/i further, the PSF 
will disappear already for a small population imbalance, 



since the particles form a conventional supcrfiuid. Here 
we do not find any phase separation. 

Finally we discuss the influence of population imbal- 
ance on the PSF in a trapped 3D cubic lattice using 
BDMFT+LDA. From our results shown in Fig. [TT] we 
conclude that here the physics is qualitatively similar to 
the 2D case. 



IV. SUMMARY 



that of the XY-ferromagnet in the corresponding system 
with repulsive interactions of equal magnitude. We have 
confirmed the stability of the PSF in a balanced Bose- 
Bose-mixture in the presence of the harmonic trap both 
in 2D and 3D. On the other hand, we found that even 
a small population imbalance can destroy the PSF. This 
novel PSF quantum phase can be detected in future ex- 
periments via the momentum distribution of pairs, which 
shows signatures of the pair condensate [17|, [34J . 



We have investigated low-temperature properties of 
Bose-Bose mixtures with attractive inter-species inter- 
action both in 2D and 3D optical lattices by means of 
BDMFT/RBDMFT. In particular, we found that the 
pair-superfluid is stable also for asymmetric hopping of 
the two species. We obtained the critical temperature 
of the PSF which we found to be of the same order as 
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